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We investigate the approach to the continuum limit of the spectrum of 
■ the Hermitian Wilson-Dirac operator in the supercritical mass region for pure 

ov 

gauge SU(2) and SU(3) backgrounds. For this we study the spectral flow of 
the Hermitian Wilson-Dirac operator in the range < m < 2. We find that 



(N 



the spectrum has a gap for < m < mi and that the spectral density at 
zero, p(0;m), is non-zero for mi < m < 2. We find that mi — > and, for 

> 

V) . m ^ 0, p(0; m) — > (exponential in the lattice spacing) as one goes to the 

continuum limit. We also compute the topological susceptibility and the size 

O - 

Q\ . distribution of the zero modes. The topological susceptibility scales well in 



the lattice spacing for both SU(2) and SU(3). The size distribution of the 
zero modes does not appear to show a peak at a physical scale. 



X 

S 1 I. INTRODUCTION 



Continuum gauge field theory works under the assumption that all fields are smooth 
functions of space-time. This assumption is certainly a valid one for quantum gauge field 
theories that respect gauge invariance: One should always be able to fix a gauge so that the 
gauge fields are smooth functions of space-time since the action that contains derivatives 
in gauge fields will not allow it otherwise. The space of smooth gauge fields typically has 
an infinite number of disconnected pieces where the number of pieces is in one to one cor- 
respondence with the set of integers 0. Every gauge field in each piece can be smoothly 
interpolated to another gauge field in the same piece but there is no smooth interpolation be- 
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tween gauge fields in different pieces. This is the case for U(l) gauge fields in two dimensions 
and SU(N) gauge fields in four dimensions. 

In lattice gauge theory, gauge fields are represented by link variables U^ipc) that are 
elements of the gauge group. Continuum derivatives are replaced by finite differences and the 
concept of smoothness of gauge fields does not apply. Any lattice gauge field configuration, 
U(i(x) = e %A ^ can be deformed to the trivial gauge field configuration by the interpolation 
C/ M (x; r) = e %rA ^ with U^{x\ 1) = U^{x) and U^x; 0) = I. Since smoothness does not hold 
on the lattice away from the continuum limit, the space of gauge fields on the lattice forms 
a simply connected space. Separation of the gauge field space into an infinite number of 
disconnected pieces can only be realized in the continuum limit. 

In this paper, we will address the following basic question: Do we see a separation of lat- 
tice gauge fields configurations into topological classes as we approach the continuum limit? 
To answer this question, we will use several ensembles of lattice gauge field configurations 
obtained from pure SU(2) and SU(3) gauge field theory. We will use a Wilson-Dirac fermion 
to probe the lattice gauge field configuration. Our motivation is the overlap formalism |2j 
for chiral gauge theories. Topological aspects of the background gauge field are properly 
realized by the chiral fermions in this formalism and therefore it provides a good framework 
to answer the above question. The hermitian Wilson-Dirac operator enters the construction 
of lattice chiral fermions in the overlap formalism and topological properties of the gauge 
fields are studied by looking at the spectral flow of the hermitian Wilson-Dirac operator 
as a function of the fermion mass. Contrary to some other approaches to investigate the 
topological properties of lattice gauge field configurations @ we do not modify the gauge 
fields, generate by some Monte Carlo procedure, in any way. 

The paper is organized as follows. We begin by explaining in Section [TT] the connection 
between the spectral flow of the hermitian Wilson-Dirac operator and the topological content 
of the background gauge field. Possible scenarios for the qualitative nature of the spectrum 
on the lattice are presented. In Section |TTT| we present numerical results on the spectral 
properties of lattice gauge field ensembles and their behavior as the continuum limit is 
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approached. Results for the topological susceptibility in pure SU(2) and SU(3) gauge theory- 
computed using the overlap definition of the topological charge are shown in Section |V[ We 
also present results on the size distribution of the zero modes of the Hermitian Wilson-Dirac 
operator. 

II. SPECTRAL FLOW, TOPOLOGY AND CONDENSATES 

The massless Dirac operator in the continuum anticommutes with 75. Therefore, the 
non-zero imaginary eigenvalues of the massless Dirac operator come in pairs, ±zA, with ij) 
and 75^ being the two eigenvectors. The zero eigenvalues of the massless Dirac operator 
are also eigenvalues of 75. These chiral zero modes are a consequence of the topology of 
the background gauge field. It is useful to consider the spectral flow of the Hermitian Dirac 
operator: 

H(m) = 75 (7A-H (1) 

The non-zero eigenvalues of the massless Dirac operator combine in pairs to give the following 
eigenvalue equation: 

H(m)x± = A±(m)x± = ±VA 2 + m 2 x± ■ (2) 

X± are linear combinations of if) and 75?/'. The eigenvalues A±(m) of these modes never cross 
the x-axis in the spectral flow of H(m) as a function of m. The zero eigenvalues, r j ll D^(j)± = 
with 750± = ±0-t result in 

H(m)</>± = =Fm0± (3) 

These modes, associated with topology, result in flow lines that cross the x-axis. A positive 
slope corresponds to negative chirality and vice- versa. The net number of lines crossing zero 
(the difference of positive and negative crossings) is the topology of the background gauge 
field. Global topology of gauge fields cause exact zero eigenvalues at m = 0. In addition, 
one can have a non-zero spectral density at zero. In an infinite volume in the continuum, the 

3 



spectrum is continuous and p(X; m)dX is the number of eigenvalues in the infinitesimal region 
dX around A. The spectral gap X g (m) defined as the lowest eigenvalue at m is equal to \m\. 
The spectral density at zero, p(0; m), can be non-zero only at m = indicating spontaneous 
chiral symmetry breaking in a theory like QCD. The continuum picture is shown in Fig. [I]. 
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FIG. 1. Continuum picture of the spectral gap and the spectral density at zero. 

To study the possible emergence of the above picture as the continuum limit of a lattice 
gauge theory picture, we need to have a lattice realization of H(m). It is important to 
note that we are interested in the spectral flow of a single Dirac fermion. With this in 
mind, we choose the hermitian Wilson-Dirac operator obtained by multiplying the standard 
Wilson-Dirac operator by 75: 

(B{U)-m C(U) \ 
HzM = . (4) 

V Ct(£7) -B(U) + mJ 

C is the naive lattice first derivative term and B is the Wilson term. We are interested in 
the spectral flow of Hx,(m) as a function of m. We note that m = 0, 2, 4, 6, 8 are the points 
where the free fermions become massless with degeneracies 1, 4, 6, 4, 1 respectively. Next we 



observe that H l (tti) can have a zero eigenvalue only if m > ||] . 

Spectral Gap Spectral density at zero 




FIG. 2. Possible scenarios of the spectral gap and the spectral density at zero on the lattice. 

We focus on the range < m < 2 and propose the following scenarios for the spectral gap 
and the spectral density at zero on the lattice and their approach to the continuum limit. 
Six different, but not completely independent scenarios are possible as shown in Fig. |2|. 

• On the lattice we have (a) and (i) with m c — > in the continuum limit. p(0;m c ) 
approaches the continuum limit with proper scaling taken into account. 

• On the lattice we have (b) and (ii) where (a) and (i) are the continuum limit. In this 
case, rrii — * and m 2 — > 0. In the limit we also get p(0; 0). 
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• On the lattice we first have (c) and (ii) going, at weaker coupling to (b) and (ii), 
where (a) and (i) are the continuum limit. The gap opens up at some m<i > mi at 
some coupling and afterwards the approach to the continuum is as in the previous 
scenario. 

• On the lattice we have (c) and (ii) where (c) and (i) are the continuum limit. In this 
case, nil ~~ > 0. However, the gap does not open up for m > in the continuum limit. 

• On the lattice we first have (c) and (iii), going to (b) and (ii) at some coupling. 
Afterwards, the approach to the continuum limit is again as in the second scenario. 

• On the lattice we have (c) and (iii) where (c) and (i) are the continuum limit. Here 
also m\ — > and p(0; m) = 0ifm>0. However, the gap does not open up for m > 0. 

We will show that numerical studies of the spectral flow on various ensembles favor the 
last scenario. Before we do that, we present a topological argument which will show that 
zero eigenvalues of Hx(m) can occur anywhere in the region < m < 8 fJ|. The spectrum 
of Hi(m) and — Hl(8 — m) are identical for an arbitrary gauge field background. Since zero 
eigenvalues can occur only for m > in H^(m), it follows that zero eigenvalues can occur 
only in the region < m < 8. It also follows that every level crossing zero from above in 
the spectral flow of Ht,(m) must be accompanied by a level crossing zero from below. In 
a single instanton background a level crossing zero from above at m + is accompanied by 
another level crossing zero from below at 2 > m_ > m + . The second crossing is due to 
one of the four doubler modes. Both m± will be functions of the size of the instanton p in 
lattice units. For p » a, m + pa and m_ pa 2. As p decreases, m + moves farther away 
from zero and m_ moves away from 2 and closer to m+. This motion as a function of p 
is smooth and for some value of p, m + = m_. The spectral flow changes smoothly as the 
configuration is changed slowly. As we move in configuration space the topological charge 
of a configuration changes. Tracing the spectral flow as a function of configurations shows 
that zero eigenvalues of H^(m) can occur anywhere in the region < m < 8. 
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III. SPECTRAL DENSITY AT ZERO 



In the previous section, we argued that H/,(m) can have zero crossings anywhere in the 
region mj < m < 2. Therefore the spectral gap is zero in this region on the lattice. This 
has direct implications for how the spectral density at zero behaves on the lattice. A careful 
study of the spectral density at zero has been performed on a variety of SU(3) pure gauge 
ensembles. This is done by computing the low lying eigenvalues of H^(m) using the Ritz 
functional ||. The low lying eigenvalues over the whole ensemble are then used to obtain 
the integral of the spectral density function, J A p(X'; m)d\' . A linear fit in A is made, and 
p(0; m) is obtained as the coefficient of the linear term. 

0.04 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




1.0 1.5 2.0 



in 

FIG. 3. p(0; m) as a function of m for various SU(3) pure gauge ensembles at gauge couplings 
P = 6.0, 5.85 and 5.7. 
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All ensembles show a peak in p(0; m) at some value of m near mi. There is a sharp rise 
to the peak from the left and a gradual fall from the peak on the right. There is a gradual 
rise again to a second peak at the location of the first set of doublers. The peak itself gets 
sharper and moves to the left as one goes toward the continuum limit. p(0;m) is non-zero 
for rrii < m < 2 in the infinite volume limit at any finite value of the lattice gauge coupling 
(see below), mi goes to zero as the lattice coupling approaches the continuum. p(0;m) 
approaches the infinite lattice volume limit from below as expected. We are fairly confident 
that we have the infinite volume limit estimate for p(0; m) at all the lattice spacings plotted 
in Fig. H 

In Fig. [| and Fig. ^j, we focus on the behavior of p(0; m) at a fixed m as one approaches 
the continuum limit. In Figure |4] we plot p(0; m) as a function of the lattice spacing measured 
in units of the square root of the string tension (the values for the string tension are taken 
from Ref. 0). In this figure p(0; m) appears to go to zero exponentially in the inverse lattice 
spacing. This is given some credence by plotting the same figure in a logarithmic scale in 
Fig. |5] where the data is shown for several values of m. For (3 = 5.7, the peak in p(0;m) is 
quite close to m = 1.2 as can be seen in Fig. resulting in a large value for p(0; 1.2). 

We remark that the p(0; m) plotted in Fig. [5] seem to favor a functional form fitting 
^ e -c/y/a £ or gach m The power of a in the exponent is a consequence of an empirical fit but 
the data presents substantial evidence for the following: p(0; m) in the supercritical mass 
region is non-zero for all finite lattice spacings. The approach to zero at zero lattice spacing 
is faster than any power of the lattice spacing. This shows that the last scenario presented 
in the previous section is favored by our numerical results. 
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FIG. 4. The approach of p(0; 1.7) to the continuum limit as a function of the lattice spacing in 
units of the string tension for (3 = 5.7, 5.85, 6.0, 6.1, 6.2 and 6.3. 
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FIG. 5. The approach of p(0;m) to the continuum limit as a function of for m = 1.2 

(diamonds) 1.5 (octagons) and 1.7 (crosses). For m = 1.2 only values at (3 < 6.0 are shown. 
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IV. TOPOLOGICAL SUSCEPTIBILITY 



In addition to studying p(0; ra), we also looked at the density of levels crossing zero in an 
infinitesimal range dm centered at m. In the continuum we expect levels crossing zero only 
at ra = 0, but on the lattice we find a finite density of levels crossing zero wherever p(0; ra) 
is non-zero on the lattice. The overlap formalism for constructing a chiral gauge theory on 
the lattice provides a natural definition of the index, /, of the associated chiral Dirac 
operator. The index is equal to half the difference of negative and positive eigenvalues of the 
hermitian Wilson-Dirac operator. A simple way to compute the index I is to compute the 
lowest eigenvalues of H^(m) at some suitably small ra before any crossings of zero occurred. 
Then m is slowly varied and the number and direction of zero crossings are tracked. The net 
number at some m t is the index of the overlap chiral Dirac operator. Since crossings occur 
for all values of ra in the range m\ < m < 2, the topological charge of a lattice gauge field 
configuration defined as the net level crossings in Hi(m) in the range [0, rat] will depend on 
m t . 

The topology of a single lattice gauge field configuration is not interesting in a field 
theoretic sense. One has to obtain an ensemble average of the topological susceptibility and 
study its dependence on ra. This has been done on a variety of ensembles and the results 
show that the topological susceptibility is essentially independent of ra in the region to the 
left of the peak in p(0; ra). A detailed study of the SU(3) ensemble at (3 = 6.0 on a 16 3 x 32 
lattice presented in Figure |^ illustrates this point. In the first line is shown the density of 
zero eigenvalues p(0; ra) and the number of crossings in each mass bin. We see that p(0; m) 
rises sharply in m, then falls to a nonzero value where there is a small number of levels 
crossing zero. In the second line of Figure || we show the size of the zero modes p z (m). We 
define a size of the eigenvector associated with the level crossing zero mode as 

Pz(m) = m^- /(t) = £tr(0t(f, t )0(£,t)) 

^ /max g 

where 4>{x,t) is the eigenvector of Hl at the crossing point and f max is the maximum value 
of f(t) over t. Another definition based on the second moment of f(t) was used in Ref. |§. 
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We should emphasize that we look only at the sizes of eigenmodes that cross, and only close 
to the crossing point. Only then can we expect to get a good estimate of the localization 
size inspired by the 't Hooft zero mode. The modes are large near mi where p(0; m) is large, 
then p z (m) drops sharply to about 1 or 2 lattice spacings and stays there up to m = 2. We 
see that the corresponding topological susceptibility rises sharply when p z {m) is large for m 
near m\ and then it is quite stable when p z (m) is small. This result shows that while the 
index, J, of the field is m dependent, the topological susceptibility, x ( a physical quantity) 
is independent of the contribution from the small modes for m > 1. 
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FIG. 6. Detailed study for = 6.0, 16 3 x 32. 
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To further clarify the relative contribution of the zero modes, in the last line of Figure |B| 
the zero mode size distribution is plotted as a function of p z . In the adjacent graph, the 
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topological susceptibility, x-, nere denned by the contribution of zero modes of size p z and 
larger, is stable when p z < 2. Hence, the small modes do not affect the estimate of \ even 
though there is an abundance of such modes. Our estimates of \ are shown in Table | where 
we use the string tension value yfcj = 440 MeV to set the scale. Our results are in rough 
agreement with other groups |J and also show good evidence for scaling. 



TABLE I. Topological susceptibility and parameters for SU(3) and SU(2). 
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V. SIZE DISTRIBUTION OF ZERO MODES 



Studies in smooth gauge field backgrounds on the lattice have shown that single instan- 
tons result in a single level crossing zero at some m in the region [0, 2] where the shape of the 
mode at the crossing is a good representation of the shape of the instanton Q. Similarly, 
there is a pair of levels crossing zero (one from above and another from below) when the 
gauge field background has an instanton and anti-instanton. Again, the shape of the modes 
at the crossing points are good representations of the instanton and anti-instanton. This 
motivates us to look at the size distribution of the zero modes of the Hermitian Wilson-Dirac 
operator on the lattice. Since lattice gauge fields generated in typical Monte Carlo simu- 
lations are rough, the correspondence between zero modes and topological objects might 
be questionable - the existence of topological objects (a collection of instantons and anti- 
instantons) "underneath" the typically large quantum fluctuations is somewhat questionable 
as well. However, levels crossing zero contribute to the global topology and an analysis of the 
size distribution of the zero modes is therefore interesting. Such size distributions are shown 
in Figure [7] for gauge group SU(2). All distributions show a sharp rise at small sizes due to 
the abundance of small zero modes that occur in the bulk of the region ni\ < m < 2. These 
modes do not affect the computation of the topological susceptibility and can be viewed as 
being due to the ultra-violet fluctuations in the gauge field background. If we eliminate the 
small modes from the distribution, the size distribution at (3 = 2.4 on the 16 4 lattice shows 
some evidence for a broad peak around p z = 0.6 fm. One should keep in mind that the box 
size is roughly 1.92 fm and the peak is occurring at a value which is roughly a third of the 
box. It is tempting to explain this peak as a finite volume effect. Some support for this 
explanation is provided by looking at the distributions at j3 = 2.5 and f3 = 2.6 on a 16 4 
lattice. These boxes are now roughly 1.38 fm and 0.98 fm, respectively. After discarding 
the small zero modes, both the distributions show a broad peak at roughly p z = 0.45 fm 
and p z = 0.3 fm, respectively. As in the (3 = 2.4 case these peaks occur at roughly a third 
of the box size and the magnitude of the peak is larger as one goes to weaker coupling for 
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a fixed lattice volume. This is quite consistent with the peak being a finite volume effect. 
In Fig. H we show all the SU(2) size distributions together plotted in lattice units. There is 
evidence for a broad peak at roughly 5 lattice units - roughly a third of the lattice box size. 
Therefore, we conclude that the size distribution of zero modes does not show evidence for a 
peak at a physical scale even after we remove the small modes which are most likely lattice 
artifacts. We have to conclude that it is not possible to relate the size analysis of the zero 
modes carried out here to a size distribution of topological objects as it is postulated for the 
instanton liquid model of QCD JlO| . 
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FIG. 7. Size distribution of zero modes on various SU(2) ensembles. 
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FIG. 8. Size distribution of zero modes on the SU(2) ensembles plotted in lattice units. The 
crosses correspond to (3 = 2.4, the diamonds (3 = 2.5, and the squares (3 = 2.6. All lattice sizes are 
16 4 . 



The conclusions remain the same for the various SU(3) ensembles that we studied. The 
size distributions are plotted in Figure [5| The (3 = 5.7, 5.85 and 6.0 ensembles come from 
lattices with linear extent roughly equal to 1.4 fm, 1.04 fm and 1.6 fm, respectively. Clearly 
the size distribution on the (3 = 5.85 ensemble suffers strongly from finite volume effects 
whereas the (3 = 6.0 ensemble is not affected as much. We should remark that we do not see 
any evidence for a finite volume effect in the computation of the topological susceptibility. 
This is probably because the size distribution of the individual zero modes is not that relevant 
for the global topology which only depends in principle on the net number of level crossings 
and not on the size and shape of these crossing modes. 
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VI. DISCUSSION 



A probe of pure lattice gauge field ensembles using Wilson fermions has revealed that 
the gauge fields are not continuum like on the lattice at gauge couplings that are typically 
considered to be weak. If they were continuum like, we should have seen evidence that 
p(0; m) is non-zero at a single value of m or in a region in m that is of the order of the 
lattice spacing. Furthermore, we should have seen a symmetry in the spectrum at values 
of m on either side of the point (or region) where p(0;m) is non-zero. Instead, we found 
that p(0; m) is non-zero in a region m x < m < 2. In the continuum limit, there is evidence 
that mi goes to zero and that p(0; m) goes to zero away from m = 0. However, the spectral 
distribution does not show evidence for a symmetry asm^ — m. 
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FIG. 9. Size distribution of zero modes on various SU(3) ensembles. 
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A remark on the approach of p(0; m) to the thermodynamic limit at a fixed lattice gauge 
coupling is in order. For a small lattice, with linear size of the order of the extent N t for 
which the finite temperature deconfinement transition occurs at the given gauge coupling, 
the number of very small eigenvalues of the Wilson-Dirac operator will be essentially zero. 
When the lattice size is increased this number will grow rapidly, leading, for a while, to 
a rapid increase in the extracted estimate of p(0; m) and then leveling off at the infinite 
volume value of p(0; m). We found that this happens for a linear size about twice the extent 
N t mentioned above. 

The density of level crossing zero modes, dN / dm, of the Hermitian Wilson-Dirac operator 
is in accordance with the behavior of p(0; m). In spite of a large number of levels crossing 
zero in the bulk of m x < m < 2, we found that the topological susceptibility is unchanged by 
these small localized modes. We therefore interpret them as due to ultra-violet fluctuations. 
The size distribution of the zero modes is dominated by these small modes. However, the 
distribution, after we remove these small modes, does not show any clear peak at a physical 
scale. Some broad peaks we see in are explained as a consequence of finite volume effects. 

We finally remark that all our studies in this paper have been on pure gauge field ensem- 
bles. There is a prediction for the spectral distribution of the Wilson-Dirac operator in full 
QCD using a continuum chiral lagrangian The prediction resembles the second scenario 



presented in section II and is qualitatively different from the result we have obtained for pure 
gauge ensembles. It would be interesting to test this prediction by numerical simulations of 
full QCD with Wilson fermions in the supercritical region. There is, however, a technical 
problem in using standard Hybrid Monte Carlo type algorithms for such simulations: the 
system will be locked in a single topological sector (with topology defined as half the dif- 
ference of negative and positive eigenvalues of the hermitian Wilson-Dirac operator at the 
supercritical mass, m^, where the simulation is carried out). This is due to the fact that a 
change in topology will require a change of net level crossings in the region < m < ma- 
However, the spectral flow has to be smooth as we update the configurations using classical 
dynamics in HMC type algorithms. Therefore, at some point in the change of topology the 
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level crossing would need to occur at m<2, but such a configuration has a vanishing fermion 
determinant and hence can not be reached. Modifications of the HMC algorithm to circum- 
vent this problems would need to be developed before a study of the spectral flow in full 
QCD in the supercritical region can be attempted. 
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